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SECTION  I 

FORMULATION  AND  SOLUTIONS 

An  approximate  theory  to  define  the  stress  field  within  a composite 
laminate  has  recently  been  developed  [Ij.  It  is  the  purpose  of  the  present 
communication  to  derive  the  solution  for  a significant  class  of  stress  con- 
centration boundary  value  problems  in  composite  laminates,  namely,  the 
free  edge  problem  [2,  3],  based  upon  this  new  theory. 

As  discussed  in  Reference  1,  we  consider  a symmetric  laminate  in  which  each 
layer  is  reinforced  by  a system  of  parallel  fibers  oriented  at  an  angle  0 with 
the  x-axis  (see  Fig.  1 of  [1])  where  the  origin  of  coordinates  is  located  at 
the  center  of  the  laminate.  The  body  is  subjected  to  forces  applied  only  on 
the  ends  such  that  a constant  axial  strain  e = e is  imposed.  Each  layer  is 
also  under  the  influence  of  expansional  strains  e^  (i  = 1,2,  3,6)  which  we 
assume  are  constants.  We  also  assume  perfect  bonding  between  adjacent 
layers.  Hence,  the  stress  field  in  this  class  of  problems  is  a function  of 
y and  z alone.  The  laminate  consists  of  N layers  which  are  identified  by  the 
index  k(k  = 1,2, •••N).  As  in  flj,  we  shall  drop  the  index  k except  when 
it  is  needed  for  clarity. 


To  begin,  we  define  the  deformation  measures  e.,  k , and  expansion- 
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where  lower  case  Latin  subscripts  have  the  range  1,2,  3,  4,  5,  6 and  Greek 
subscripts  assume  the  values  1,2,  3,6.  Letting  S_  represent  the  (mono- 
clinic)  compliance  matrix  as  given  by  eq.  (12)  of  [1],  we  make  the  further 
definitions 


A. 

ij 


B.. 


[S..  + \ 6,.  6,. 
ij  5 3i  3j  SS-* 


[S..  + I 6,.  6 s 

^ ij  7 3i  3j  33-' 


(2) 


where  6.  is  the  Kronecker  delta  and  the  superscript  -1  stands  for  matrix 
inverse.  Throughout  this  work,  summation  over  the  range  of  repeated  sub- 
scripts, but  not  superscripts,  is  implied.  We  may  now  invert  the  consti- 
tutive relations,  (25)  of  [1),  to  get 
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where  standard  contracted  notation  has  been  employed  in  the  first  two 
expressions,  following  the  subscript  relations  given  in  (8)  and  (9)  of  [1]. 

Since  the  stresses  in  the  present  class  of  boundary  value  problems 
only  depend  on  y and  z,  it  follows  that  the  force  and  moment  resultants  and 
the  interlaminar  stresses  will  only  depend  on  y.  Hence,  from  (3),  the  de- 
formation measures  are  functions  of  y alone  and  by  use  of  eqs.  (1),  it  can 
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be  shown  that  the  most  general  form  of  the  weighted  displacements  is  given 
by 

u = U(y)  + CjXy  + c^x 
2 

- 

V = V(y)  - — — + c^x 

_ 2 
hw  = W(y)  - 6c^xy  - 3c^x  + 3c^x 

(4) 

U*  = l/f(y)  + c^x 
V*  = fi(y)  + c^x 
w*  r 0(y) 

A 2 

hw  = y(y)  - 2c  xy  - c X + c x 
5 6 4 


where  U X are  arbitrary  functions  of  y and  c c,  are  constants. 

1 6 

We  should  recall  that  eqs.  (1)  - (4)  must  be  written  for  each  layer. 


We  now  substitute  eqs.  (4)  into  (1),  then  into  (3),  and  finally  into 
the  equilibrium  equations,  (26)  of  [ 1 ] , in  which  the  x dependence  is  dropped 
to  establish  the  following  relations 
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which  are  also  valid  within  each  layer.  The  remaining  field  equations,  the 
interface  continuity  conditions,  are  given  by  substituting  (4)  into  (27)  of  [ IJ, 
which  gives 
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(k)  (k+n  j (k  + 1)  (k)  2 
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and 


(k  + 1)  (k) 


(k+1) 
s , 


I 


(k) 


(k  + 1) 

P, 


(k) 
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(7) 


for  k = 1,2  • • • N-  1 . 


Due  to  the  symmetric  lamination  geometry,  the  interlaminar  shear 
stress  components  and  the  z displacement  component,  w,  all  vanish  on  the 
central  plane  z = 0.  We  shall  take  advantage  of  these  conditions  by  consid- 
ering only  the  upper  half  of  the  laminate,  i.e.,  z5?0.  Incorporating  the 
traction- free  conditions  on  the  upper  surface,  our  boundary  conditions  on 
the  upper  and  lower  surfaces  become 
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Since  eqs.  (6)  and  (8)  must  be  satisfied  for  all  values  of  x,  it  follows  that 
(k)  (k)  (k) 


k = 1,2 
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and 
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so  that  Cj,  c^  and  c^  are  the  same  for  every  layer. 
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We  now  turn  our  attention  to  the  edge  boundary  conditions,  which 

require  consider. .tion  ofN.N  ,V,M,M  ,s,  and  s_  for  each  layer 

y xy  y y xy  1 Z 

on  y = ± b,  since,  as  discussed  in  [1],  only  edge  tractions  are  imposed  in 
the  present  class  of  boundary  value  problems.  However,  all  these  functions 
cannot  be  independently  prescribed  because  of  the  consequences  of  inter- 
face continuity  and  global  equilibrium  of  the  entire  laminate.  That  is,  the 

interface  continuity  conditions  given  by  the  second  of  (7)  prohibit  arbitrarily 

(k)  (k)  (1)  (Nf 

prescribed  values  of  s and  s.  . Furthermore,  s and  s have  already  been 

lb  lb 

specified  by  (9)  for  all  values  of  y.  These  relations,  in  conjunction  with  the 
second  equilibrium  equation,  see  (26)  of  [Ij,  can  be  used  to  establish  the 
result 


which  requires  that 


(k) 
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y.  y 
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(12) 


N (k)  N (k) 

E E = 0 

k = 1 y k = 1 ^ 

(k) 

Therefore,  only  2N-  1 values  of  N can  be  arbitrarily  prescribed  on  the 

y (k) 

edges  y = ± b.  We  can  make  the  same  statement  regarding  N since  an 
equation  of  the  form  (13)  can  be  derived  in  similar  fashion  for  this  function. 
Hence,  the  edge  boundary  conditions  may  be  expressed  as 
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= N (- 
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b)  = 3^{±h) 

= 0 

(k  = 1, 2-  . . N-  1 ) 

(14) 

The  present  boundary  value  problem  therefore,  consists  of  differential 
equations  (5)  and  (6)  subject  to  the  boundary  conditions  (8),  (9),  and  (14). 
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The  general  solution  for  each  dependent  variable  consists  of  the  sum 
of  two  parts;  i)  a complementary  solution  defined  by  the  homogeneous  form 
of  (5)  - (9)  and,  ii)  a particular  solution.  In  the  particular  solution  (denoted 
by  subscript  P) , the  only  nonvanishing  functions  are  given  by 
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a . 
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where  a.  (i  = 1,2)  are  constants  given  by  substituting  (15)  into  (5),  (6)  and 
i 

(8)  to  get 
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where  we  have  put 


2c, 


^1  ^3 


17) 


and  e is  the  applied  axial  strain.  The  second  condition,  c^  = 0,  follows 
from  the  procedure  leading  to  ( 1 5)  and  (16),  while  c^  is  set  to  zero  since  it 
has  no  effect  on  the  stress  field  here. 


Since  the  field  equations  are  linear  differential  equations  with  constant 
coefficients,  the  complementary  solution  (subscript  H)  for  each  dependent 
variable  consists  of  a series  of  terms  of  the  form 
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•4l 


F e 
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<k)  (k) 

where  f stands  for  any  of  the  dependent  variables  and  F are  constants. 

Substituting  (18)  into  the  homogeneous  form  of  eqs.  (5)  - (9)  leads  to  a 

system  of  1 3N  linear  algebraic  equations.  The  values  of  X are  determined 

by  setting  the  determinant  of  the  coefficients  equal  to  zero.  Algebraic 

expressions  for  the  expansion  of  the  determinant  were  not  written  owing  to 

the  complexity  of  such  expressions,  even  in  the  simplest  cases.  Rather, 

computer  calculated  values  of  the  determinant  for  specific  values  of  A,  as 

discussed  later,  were  employed  to  continue  the  analysis.  Hence,  it  is  not 

possible  to  exhibit  the  mathematical  details  of  the  remaining  steps  in  the 

solution.  Therefore,  this  phase  of  the  work  will  be  presented  in  a 

descriptive  fashion. 

Although  our  determinant  is  too  complex  to  evaluate  in  general  terms, 
we  can  examine  the  nature  of  its  polynomial  (in  X)  expansion  in  some  detail 
for  fairly  small  values  of  N and  extend  the  results  by  induction*.  We  can 
also  develop  a procedure  to  search  for  terms  in  the  determinant  in  which 
the  highest  and  lowest  powers  of  X occur.  Proceeding  in  this  fashion  for 
values  of  N = 1 (in  which  case  the  interface  continuity  conditions  are  dropped) 
2,  3,  and  4,  we  make  the  following  observations;  i)  Only  even  powers  of  X 

4 

occur  in  the  determinant,  ii)  The  lowest  power  of  X is  A , and  iii)  The 
highest  power  of  X is  A^^^  Although  these  results  are  not  perfectly 
rigorous,  no  contradictions  have  been  found.  Furthermore,  the  number  of 
X roots  is  consistent  with  the  number  of  edge  boundary  conditions,  as  dis- 
cussed later. 

An  exception  to  the  form  (18)  occurs  with  the  appearance  of  repeated 
roots  for  X.  Since  repeated  zero  roots  always  occur,  it  is  necessary  to 
treat  the  corresponding  part  of  the  solution  separately.  In  this  part  of  the 
solution,  (18)  is  replaced  by  a third  degree  polynomial  in  y.  Representing 
the  functions  corresponding  to  the  repeated  zero  roots  by  subscripts  Ho.  we 

Although  we  have  not  done  it  here,  it  is  possible  to  develop  a program  to 
define  the  powers  of  X occurring  in  each  term  of  the  determinant  for  arbi- 
trary values  of  N by  use  of  a computer. 
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find  from  the  homogeneous  form  of  (5)-(9)  that  the  only 
functions  are  given  by 
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, k = 1,2  • • • N 
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The  constants  A and  C define  rigid  body  translation  of  the  laminate  as  a 
o o 


unit.  The  remaining  constants  in  (19)  can  all  be  expressed  in  terms  of  A 


and  C^.  Hence,  two  constants  which  effect  the  stress  distribution  have  been 
introduced  in  the  repeated  zero  part  of  the  homogeneous  solution. 


The  remaining  portion  of  the  complementary  solution  consists  of 
functions  of  the  form  (18)  corresponding  to  the  12N-  6 nonzero  values  of  A 
(we  are  assuming  that  these  roots  are  all  distinct).  These  roots,  which  occur 
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in  combinations  of  the  form  ± (a±ib),  where  b may  vanish,  were  determined 

2 

by  computing  the  value  of  the  determinant  for  6N-  2 values  of  A and  using 
these  results  to  solve  for  the  coefficients  of  the  equivalent  polynomial.  Once 
the  polynomial  has  been  established,  its  roots  can  be  determined  through 
standard  computer  routines.  In  the  usual  manner,  one  equation  in  the  homo- 
geneous version  of  (5)-  (9)  may  be  dropped  and  the  reduced  system  used  to 
relate  all  but  one  of  the  arbitrary  constants  to  the  remaining  constant  for 
each  value  of  A.  This  procedure  leads  to  12N-6  arbitrary  constants.  These 
constants,  augmented  by  and  C^,  are  evaluated  via  the  12N-4  edge 
boundary  conditions  (14).  Once  these  constants  are  defined,  the  complete 
solution  for  any  of  the  1 3N  functions  appearing  in  (5)-(9)  is  given  by 

(k)  12^6  (k)  X (k)  (k) 

f = y F e + f + f (21) 

m Ho  P 

m = 1 


where  the  last  two  terms  are  defined  by  (15),  (16),  (19)  and  (20).  The  force 
and  moment  resultants  are  now  given  by  (3)  upon  using  the  results  of  (4) 
and  ( 1) . 


A special  case  arises  in  the  solution  that  requires  separate  treatment. 

This  case  is  defined  by  the  vanishing  of  compliance  coefficients  S ,,  S_,, 

16  do 

and  well  as  expansional  strain  e^,  in  every  layer.  The  situation 

arises  when  each  layer  is  isotropic  and/or  oriented  at  an  angle  of  0°  or  90°. 

(k)  (k)  (k)  (k)  (k)  (k) 

This  leads  to  the  vanishing  of  U,  <p,  t, , t_,  N , and  M . Consequently, 

1 2 xy  xy 

dropping  the  appropriate  field  equations  and  boundary  conditions  we  find 
that  setting  the  determinant  to  zero  leads  to  only  two  zero  roots  for  A and 
8N-4  nonzero  roots.  The  number  of  boundary  conditions  (14)  reduces  to 
8N-3.  Otherwise,  the  treatment  presented  here  remains  unchanged. 


The  occurrence  of  very  large  magnitudes  of  A for  large  N in  the  present 
formulation  leads  to  values  of  e*^^  which  exceed  computer  limits.  This  in 
turn  restricts  the  values  of  N which  can  be  treated.  For  example,  N=  6 was 
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the  largest  value  permissible  for  the  properties  assumed  in  [ 1 ] . Resolution 
of  this  difficulty  is  now  being  considered. 

Numerical  results  of  the  present  solution  are  presented  in  flj. 
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